Magnetism in BN Nanotubes Induced by Carbon Substitution 
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We performed ab initio calculation on the pristine and carbon-doped (5,5) and (9,0) BN nanotubes. 
It was found that Carbon substitution for either boron or nitrogen in BN nanotubes can induce 
spontaneous magnetization. Calculations based on density functional theory with the local-spin- 
density-approximation on the electronic band structure revealed a spin polarized, dispersionless band 
near the Fermi energy. The magnetization can be attributed to the carbon 2p electron. Compared 
to other theoretical models of light-element or metal-free magnetic materials, the Carbon-doped BN 
nanotubes are more experimentally accessible and can be potentially useful. 

PACS numbers: 75.75.+a,71.15.Mb 

The discovery of room-temperature weak ferromagnetism in all-carbon system consisting of polymerized Ceo Q 
stimulated wide interest in magnetism of all-carbon materials. One year later after the above discovery, Esquinazi et 
al.fy detected a ferromagnetic signal from an oriented graphite which behaved quite differently from known magnetic 
impurities, suggesting an intrinsic origin of magnetism in graphite. In order to understand the unexpected magnetism 
in all-carbon system, some density functional theory (DFT) calculations have been carried out. Lehtinen et 
performed ab initio local-sphi-density approximation (LSDA) calculations to study the properties of a carbon adatom 
on a graphite sheet and found that this defect has a magnetic moment of up to 0.5 [1b- This was important in 
understanding the magnetism observed experimentally in graphite sheet. More recently, Ma et aL0 studied the 
magnetic properties of vacancies in graphite and carbon nanotubes. For graphite, the vacancy is spin-polarized with a 
magnetic moment of 1.0 \xb- The vacancy in carbon nanotubes can also induce magnetism, depending on the chirality 
of the nanotubes and the structural configuration with respective to the tube axis. These calculations provided some 
understanding of the observed magnetism in pure carbon systems. 

The weak ferromagnetism in pure carbon system also stimulated interest in searching for light-elements or metal- 
free magnetic materials in view of their potential applications as high temperature magnets since metal magnets lose 
their ferromagnetism at high temperatures. One of such metal- free magnets was nanographite or graphite-ribbon. 
Fujita et al. [5j performed tight binding band structure calculations on graphite ribbons with armchair and zigzag 
edges, respectively. Ribbons with zigzag edges show a sharp peak in density of states at Fermi level, indicating a 
possibility of spontaneous magnetization. Mono-hydrogenation such as zigzag-edged graphite ribbon could create a 
ferromagnetic spin structure on their edges. Based on results of first principle pseudopotential calculation, Kusakabe 
and Maruyama 6| further predicted that the zigzag edged graphite ribbon could also have finite magnetization at edges 
by hydrogenating each carbon with two hydrogen atoms while hydrogenating each carbon with a single hydrogen atom 
at the other edge. However, such a hydrogenation of carbon edges is only of interest for theorists. Another recently 
proposed model by Choi and coworkers^ is the heterostructured C-BN nanotube. They calculated the electronic 
structure of the (9,0) Ci(BN)i and C2(BN) 2 nanotubes using density functional theory and found the occurrence of 
magnetism at the zigzag boundary connecting carbon and boron nitride segments of tubes. However, in view of the 
growth conditions of carbon and boron nitride nanotubes, fabrication of the heterostructured C-BN nanotubes may 
be impractical. 

In this letter, we report results of our DFT-LSDA studies on the electronic structures of the pristine and carbon- 
doped (5,5) and (9,0) BN nanotubes. We find that the carbon substitution for either boron or nitrogen atom induces 
magnetization of the doped system, with a total magnetic moment of 1 /xg. Compared to previous models of metal- free 
materials, our proposed C-dopcd BN nanotubes are experimentally accessible as far as the structural configuration is 
concerned. 

Our calculation is base on the density functional theory and the local spin density functional formulismQ- The 
calculation was done using the SIESTA code @. The valence electrons are described by linear combination of numerical 
atomic-orbital basis set and the atomic core by norm-conserving pseudopotentials. The pseudopotentials generated 
using the Troullier and Martins scheme^(| are used to describe the interaction of valence electrons with atomic 
core and their nonlocal components are expressed in the fully separable form of Kleiman and Bylander 0, fl2j |. 
The generalized gradient approximation correction in the form of Perdew, Burke and Ernzerhof|l.lj is adopted for 
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the exchange-correlation potential. The doublc-zcta plus polarization atomic orbital basis set is employed in the 
calculation. The Hamiltonian matrix elements are calculated by charge density projection on a real space grid with 
an equivalent plane wave cutoff energy of 70 Ry. Periodical boundary condition is employed along the nanotube axis, 
and a vacuum region of at least 10 A is assumed between the nanotube and its images, to avoid tube-tube interaction. 
Special k-mesh is chosen along the tube axis according to the Monkhorst-Pack scheme01- The conjugate gradient 
algorithm 0] was adopted to fully relax the structure of the nanotube until the maximum force on a single atom is 
within 0.02 eV/A. For the calculation of carbon doped (5,5) nanotube, we use a supercell which consists of 1 x 1 x 2 
primitive cells and contains 80 atoms to avoid possible coupling between carbon images. 

The calculated band structures of pristine (5,5) and (9,0) BN nanotubes are shown in Fig. 1. The band gaps of 
4.47 eV for (5,5) and 3.81 eV for (9,0) nanotubes agree well with those of ab initio calculation by Xiang et aL[l|| The 
(5,5) BN nanotube has an indirect band gap while the (9,0) has a direct band gap. The calculation also reproduces 
the free-electron-like lowest conduct bands as found by Blase et al. These results confirm the validity of the 
application of atomic orbital basis set. We also performed spin-polarized calculation on these two pristine nanotubes 
and no spontaneous magnetization was found. 

However, when either boron or nitrogen is substituted by carbon, a spontaneous magnetization is induced. The 
band structure of the (5,5) BN nanotube with a boron atom substituted by carbon atom is shown in Fig. 2(a). We 
can see that the carbon substitution causes a lift up of the Fermi level. All bands with energies lower than -6.0 eV are 
spin degenerate and fully occupied, and thus do not contribute to spin-polarization. However, close to the Fermi level, 
a very flat band, almost dispersionless in the whole first Brillouin zone, is split. Only the spin-up branch is occupied 
and the spin-down branch is left unfilled, leading to a strong spontaneous magnetization in the (5,5) nanotube. The 
splitting between the occupied and unfilled branch is as large as 1.5 eV, The flatness of the spin-polarized band 
indicates that the corresponding electron state is heavily localized, suggesting that the magnetic moment is localized 
at the substitution site, possibly from the carbon 2p electron. Wc also investigated single carbon substitution for 
nitrogen on the (5,5) nanotube and found that it also produces a spin polarized band around the Fermi level. However, 
the Fermi level is pushed down. The net magnetic moment is about 1 [is in both substitutions. 

The band structure of the (9,0) nanotube with a nitrogen atom substituted by carbon atom is shown in Fig. 2(b). In 
this case, the Fermi level is pushed down, similar to the case of carbon substitution for nitrogen in the (5,5) nanotube. 
A spin-polarized band occurs around the Fermi level with the spin-up branch occupied and the spin-down branch 
unfilled, resulting in a total magnetic moment of 1 [Ib- Similar calculation was carried out for the case of carbon 
substitution for boron atom and it was found that the Fermi level increases and the band below the Fermi level is 
partially occupied and spin-polarized, leading to a net magnetic moment of 1 \ib- 

The spin polarization can also be clearly seen from the spin density of states (DOS), as shown in Fig. 3. For carbon 
substituted (5,5) nanotube (Fig. 3 (a)), a separate occupied majority spin peak emerges at the energy of —4.0 eV. 
Below this peak, the majority and minority spin density of states are essentially identical. For the carbon substituted 
(9,0) nanotube (Fig. 3 (b)), we also can see that near the Fermi level, the majority spin DOS exceeds that of the 
minority. However, the highest occupied spin-up branch merges with the lower valence bands in this case. 

Boron nitride has important advantages over carbon nanotube. It is far more resistant to oxidation than carbon 
and therefore more suitable for high-temperature applications in which carbon nanostructures would burn. Unlike 
carbon nanotubes, BN nanotubes are insulators, with predictable electronic properties that are independent of their 
chiralities. In addition, magnetic nanostructures are of scientifically interesting and technologically important, with 
many present and future applications in permanent magnetism, magnetic recording and spintronics. The carbon 
substituted BN nanotubes, with conduction electrons that are 100% spin polarized due to the gap at the Fermi level 
in one spin channel, and a finite density of states for the other spin channel, can be an ideal half metallic material and 
can be useful for spintronic applications, such as tunneling magnctorcsistancc and giant magnctoresistance elements. 
Upon geometry optimization, all doped nanotubes were spontaneously spin polarized and no visible deformation or 
buckling in their structures can be observed. Contrast to the models proposed by Kusakabe and Choi Q in 
which the spin polarization only occurs for certain particular structural configurations, the spin polarization in BN 
nanotubes is induced by substitution with arbitrary atomic configuration, indicating easy accessibility of experimental 
synthesis. 

In conclusion, we have performed first principle pseudopotential calculation with double-zeta plus polarization 
atomic-orbital basis set to study the effects of carbon substitution on BN nanotubes. Our results show that single 
carbon atom substitution for any atom results in a polarized fiat band and a sharp peak in the majority density of 
state below the Fermi level. The spontaneous spin-polarization is independent of site of substitution and chirality of 
the nanotubes. Compared to other metal-free magnets previously proposed, the carbon induced magnetization in BN 
nanotubes is experimentally accessible and the system can be potentially very useful. 
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FIG. 1: The band structure of the pristine (5,5) (a) and (9,0) (b) BN nanotubes. 
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FIG. 2: (a) The band structure of (5,5) nanotube with a boron atom substituted by carbon, (b) The band structure of (9,0) 
nanotube with a nitrogen substituted by carbon. Solid lines and + represent the bands for spin-up and spin-down electrons 
respectively. The Fermi level is denoted by the dotted line. 
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FIG. 3: Majority and minority spin densities of states of the (5,5) (a) and (9,0) BN nanotubes. 



